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Abstract. A system-size expansion method is incorporated into the Doi-Peliti 
formalism for stochastic chemical kinetics. The basic idea of the incorporation 
is to introduce a new decomposition of unity associated with a so-called Cole- 
Hopf transformation. This approach elucidates a relationship between two different 
Langevin equations; one is associated with a coherent-state path-integral expression 
and the other describes density fluctuations. A simple reaction scheme X <± X + X is 
investigated as an illustrative example. 

PACS numbers: 05.10.Gg,05.40.-a,82.20.-w 
1. Introduction 

System-size expansion methods were developed in the analysis of master equations 
for stochastic chemical kinetics [1-3]. In these methods, under the assumption of 
an extensive property of macroscopic quantities, fluctuation effects are systematically 
taken into account by assuming the inverse of a dimensionless system size to be a 
small perturbation parameter. The leading-order description of a system is provided 
as a deterministic equation for a density variable due to the law of large numbers, and 
macroscopic fluctuations are described by a Langevin equation which reduces to the 
deterministic equation in the limit of an infinite system size. 
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The Doi-Peliti formalism is known as another standard technique, where master 
equations are equivalently expressed in terms of a set of creation and annihilation 
operators [4,5]. Since the original proposal, many models have been analyzed within this 
formalism, as seen in a list of papers quoted in Ref. [6]. In particular, a coherent-state 
path-integral expression is useful for theoretical study, because it takes a form similar 
to that for a continuum field theory [6-8]. 

Then, in some cases including an example discussed in this paper, the dynamical 
action of a coherent-state path-integral expression becomes equivalent to that of a 
Langevin equation. However, since this Langevin equation does not coincide with that 
obtained by system-size expansion methods, it does not directly describe the stochastic 
time evolution of a density. Note that this difference is not recognized well; indeed, 
as presented in a review paper [9], a Langevin equation obtained by a coherent-state 
path-integral expression [10] has been analyzed in the problems of fluctuating front 
propagation of some reaction-diffusion systems as if it describes the time evolution of a 
density. (See however a different argument in Ref. [11].) 

In order to resolve this confusing situation, in the present paper, we develop a 
system-size expansion method for a coherent-state path-integral expression in the Doi- 
Peliti formalism. Our consideration consists of three steps. In the first place, we 
express statistical quantities in terms of a variable that is straightforwardly introduced 
in the coherent-state path-integral expression. In addition to the presentation of such 
an expression, which was already discussed in Refs. [8, 10], we shall provide a clear 
argument on the indirect correspondence between density fluctuations and trajectories 
of the Langevin equation associated with the coherent-state path-integral expression. 
At the second step, we seek for a simple expression of density fluctuations in terms 
of coherent states. Although it has been argued that such an expression would be 
realized by a so-called Cole-Hopf transformation [12-15], the description involves non- 
transparent procedures, as will be indicated in the beginning of section HI Our main idea 
in the present paper is to introduce a new decomposition of unity associated with the 
Cole-Hopf transformation, by which we obtain a new path-integral expression. Finally, 
this expression enables us to perform a system-size expansion, and it directly leads to 
the Langevin equation describing density fluctuations. 

This paper is organized as follows. In section [2J we introduce a model we study 
within the Doi-Peliti formalism. In order to make our presentation as instructive and 
transparent as possible, we focus on a simple reaction scheme X *=± X + X in a single 
box under an assumption that diffusion processes are sufficiently fast. In section [31 
we review a coherent-state path-integral expression, and derive the first type Langevin 
equation. In section HJ we propose a new decomposition of unity, by which we derive 
the second type Langevin equation. In section [5l we briefly discuss effects of diffusion. 
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2. Model 

Let Q be a volume of a box in which particles are confined. We consider a chemical 
reaction X + X ^ X; the duplication reaction X — > X + X occurs with ratio a for each 
particle, and its backward reaction X + X — > X occurs with ratio /? when two particles 
are close to each other less than a reaction length tq. We assume that a diffusion time 
is much smaller than typical time scales of the reaction: a -1 and and ignore the 
effect of diffusion for the time being (we will briefly discuss it in section [5]) . Then, 
since the probability of finding another particle within a reaction region for a given 
particle is 47rrQ/(3f2), the rate of the reaction X + X — > X for each pair of particles 
becomes PAirr^/ (30) . In the argument below, a dimensionless quantity Anr^/fifl) is 
replaced with for notational simplicity. Now, let P n (t) be a probability of finding 
n particles at time t. The time evolution of P n (t) is described by the master equation 

j t P n (t) = a[(n - l)P n _i(*) - nP n (t)} + £[n(n + l)P„+i(t) - n(n - l)P n (t)] , (1) 
where P-i(t) = 0. 

State of the system is specified by a vector |n), which represents a situation that 
there are n particles in the system. Let V be a vector space spanned by an orthogonal 
set: V = { |n) | n = 0, • • • , oo}. A remarkable idea of the Doi-Peliti formalism [4, 5] is 
that a series of infinite number of P n (t) is collectively treated as a single vector in V: 

oo 

i^)> = E p »wi n >- ( 2 ) 

n=0 

Then, the master equation (TjQ) is equivalently rewritten in a compact form 

J^(f)> = -^I^(O), (3) 

where is an infinite-dimensional matrix. This matrix is expressed in terms of two 
matrices a and a) defined by a\n + 1} = (n+ 1) |ri) and a'\n) — \n+ 1) for any n (with 
o|0) = 0): 

- a (a f - l)a f a- ^(1 -a^o+o 2 . (4) 

The two matrices satisfy the bosonic commutation relations [a, a'] = 1 and [a, a] = 
[a^a^] = 0. Since and a are similar to the creation and annihilation operators, and 
H is an analog of the Hamiltonian, the Doi-Peliti formalism is also called the second- 
quantization method. 

For any \v\) and \v2) in V, the inner product (vi\v2) is naturally defined by noting 
(m\n) = n\ 5 TO>n . There is a special vector \V) such that (n\V) = 1 for any n, which 
is employed to express the expectation value of observables. This vector is explicitly 
written as \V) = e at |0). Let A{x) be a polynomial of a variable x. We consider the 
expectation value of A(n) at time t — r, which is denoted by (A(n(r))). By using 
the vector |"P), we write (A(n(r))) = J2^=o A{n)P n {r) = (V\A(a) 'a)\tp{j)) . In order to 
simplify this expression further, we define a polynomial A(x) associated with A(x) by 
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the relation (V\A(a^a) = (V\A(a). That is, A(a) is obtained by the normal ordering 
of A{a*a). As an example, we demonstrate for the simplest non-trivial case A(x) = x 2 . 
Since {V\(a)a) 2 = (V\aaJa = (V\(a 2 + a), one finds A(x) = x 2 + x. More generally, by 
mathematical induction, one can prove that A(x) = x k for A(x) = x(x — 1) • • • (x — k + 1) 
with any integer k. The final expression of the expectation value of A(n) at t — r 
becomes 

(A(n(r))) = (V\A(a)\iP(r)} . (5) 
3. Coherent-state path-integral expression 

With the "second-quantized" expression of the master equation ([3]), it is quite natural to 
consider the path- integral form of the expectation value ([5]). In this section, we explain 
a procedure to obtain a path-integral expression by using coherent states within the 
model (TjJ. The result was already obtained in Ref. [10], but we shall provide a slightly 
different presentation of the derivation. 

Coherent states are eigenstates of the matrix a. They are in general written as 

|^)=e 2at |0) (6) 

with a complex number parameter z being the eigenvalue: a\z) = z\z). The 
corresponding bra vector is denoted by (z\, that is, (z\ = (0|e 2 * a , where z* is the 
complex conjugate of z. The following decomposition of unity plays an essential role in 
constructing a path-integral expression with the coherent states: 

oo oo oo „ _. 2 

where we have used 

6 n , m = I ^e-' 2 'V"V\ (8) 
J irnl 

with the integration measure d 2 z = d(Re z)d(Im z). See Refs. [6] and [8] for the path- 
integral expression on the basis of (J7J). 

Here, we notice a different decomposition of unity 

1= f°d0 r^|0>(-Me-^, (9) 

where we have used the formula 

n\ J-oo V <W n! J.oo 2tt 

Note that and ip are real variables. The path-integral expression on the basis of 
the decomposition of unity (Q might be less familiar, but we remind readers that this 
expression was presented in the seminal paper by Peliti [7]. It should be equivalent to 
the standard one obtained by using (j7j). Indeed, when we employ the expression with 
(4>, iip), we can always move to the standard expression by formally replacing (0, 'up) with 
(z,z*). In addition, as far as we understand, the expression obtained by flU]) is more 
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transparent in the sense that the integration can be considered in an explicit manner. 
From these reasons, we review the coherent-state path- integral expression by using Q. 
The time-integration of ([3]) yields 

|V(r)> = hm o [eM-HAt)Y^\m) • (11) 

We assume that an initial condition obeys a Poisson distribution with its average n . 
The initial state \ip(0)) is then expressed as 



Pois) = ]Te-"°^|z) = e" ^- 1 ^) . 



i=0 



For the moment, we focus on the case 



ano 



Uq > . 



(12) 



(13) 



and the other case will be discussed at the end of this section. By inserting the identity 
(ED into (ED with (ED and pj, we express (A(n(r))) as 



(A(n(r))) 



lim 

Ai^O 



n 

A=0 



dipt 



2tt 



(V\A{a)\<fH) 



x 



;Q(-i^|e-^ A U-A/ 



)e- i0t 



•ft 



.t=At 



-iip \n ; Pois) 



(14) 



where the time index t runs from zero to r in steps of At. (We assume r = NAt with 
integer N.) Since each term of H in (jlj) is arranged in normal order, we derive 

(-i^|e-^ A *|0i_ At ) = (-i^t-A*) exp(-#(i^, &-At)At) , (15) 

where H(i<p t , 4>t-At) is defined by simply replacing a) and d in H with fy? t and <fit-At, 
respectively. We also have (-iipt\<j>t-At) = exp(iip t 4>t-At) , (V\A{a)\4> T ) = (l\(j) T )A((j) T ) = 
A((j) T )e^ T , and (— i<p \n ; Pois) = exp (iy^o — ^o) • By substituting these results into 
([141) . we obtain a path-integral expression 



(A(n(r))) 



lim 

At-»0 



n 

v t=0 



dipt 



2tt 



^(0 T )exp[-S({i^},{0})] , 



(16) 



where the "action" S is calculated as 
^({i^}, {0}) = + iVo(0o - n ) + n 



t=At 



h-At 



At 



(17) 



Although S is a complex valued functional as it is, the integration with respect to (ft in 
(fT6l) yields a real valued functional of (p. Here, the integration of (ft in (flBT) is carried 
out from — oo to +oo along the real axis. We assume that the integration path can be 
shifted to the straight line from — oo — i to +oo — i with keeping the integration value 
unchanged. This shift is equivalent to the replacement of ip t with (p t — i, where (p t is a 
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real variable. Then, the action (TTT|) becomes 
^({i^}, {<£}) = i<^o(0o - n ) 



t=At 



m- 



H-At 



(3 



At 



- " ((m) 2 + m) <k-M + ^ (m + (m) 2 ) 4> 



2 

t-At 



(18) 



Next, we linearize the terms quadratic in (p t by introducing an auxiliary real variable y t : 
exp [(i^) 2 (afo-M ~ Wm 2 t-At) (At)] 



cxp 



--y\ + Mp t Jl (a4> t -At ~ GW#-At) (At) 



1/2 



Vt 



(19) 



Note that this identity is valid only when 

«^-A*-(/3/fi)0 2 _ At >0. (20) 

We conjecture that trajectories that do not satisfy this condition do not contribute 
the path- integration in ffT6]) under the condition fflBl . See a discussion below f f24|) . 
Finally, by employing the Fourier transformation formula of Dirac's delta function, we 
can rewrite (fTBl) as 



<^Wr))>=Um/n / 

\ +— A + J — 



t=At J -°° 



dyt 
. ( 

V2^ 



A{(j) T ) 



Y[6m<t> t ,<f) t _ At ,y t )) 



t=At 



,(21) 



<f>0=no 



where 



®(<f>t, 4>t-At, Vt) =<Pt~ (f>t-At - ( 0«f)t-At ~ o 



f 3 ^ 



(At) 







t-At 



- \ l 2 ( a&t-At 
The equation $ = in the limit At — > is expressed as 

■V, 



(Aty"y t 



^■o no - ^o 2 + ^i2 I n< 



(22) 



(23) 



where the quantity r)(t) defined as the continuum limit of y t /v At is a white Gaussian 
noise that satisfies (r](t)r](t')) = 5(t — t'), and the symbol ■ in front of 7] represents the 
product in the sense of Ito. Therefore, the formula ( !2~Tl) with ( l22l) is written as 

(A(n(r))) = (A{<Kt)))+ , (24) 

where ( )^ represents the expectation value over noise sequences of the Langevin 
equation ( l23l . It should be noted that trajectories (pit) generated by ( l23l) with the 
initial condition 0(0) = no satisfy acf) — (/3/Q)<p 2 > for any time t as discussed in 
Ref. [10]. Then, one may obtain ( 1181) starting from ( 1231) . This implies that the conjecture 
mentioned below (1201) is valid. 

We choose A(n) = n/Q in (124")) as the simplest example. By setting u(t) = <p(t)/Q, 
we have 

'n(r) 



(«(t)>„ 



(25) 
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where ( ) represents the expectation value over noise sequences of the Langevin equation 



d _ 2 2(au — (3u 2 ) , n „. 
-u = au-(3u 2 + ^ -JL-L. V . (26) 

Surprisingly enough, the average behavior of the density n/Q. is exactly described by 
the Langevin equation (I26p for any Q. The Langevin equation (j26p supplemented with 
a diffusion term is called the stochastic Fisher Kolmogorov-Petrovsky-Piscounov (s- 
FKPP) equation, which appears in several research fields [9] (for example, see Ref. [16] 
for applications to the problems of high-energy hadron scattering). Historically, the 
s-FKPP equation for the chemical reaction X ^ X + X was derived on the basis of the 
coherent-state path-integral expression with the complex parameterization (z, z*) [10] 
in studying non-trivial noise effects for front-propagation [17,18]. Here, it is worthwhile 
noting that does not correspond to the number of particles, as we can see the difference 
between the two polynomials A and A that appear in (f24j) . It is only when we take 
the simplest case A(x) = x discussed in ( 1251) that A(x) becomes equivalent to A(x). 
Therefore, one must not interpret that (1261) is the equation for a fluctuating density. 
Indeed, ( )26l) is different from that obtained by system-size expansion methods for the 
master equation ([T]). (Compare fl2"B"l) with fl5T)|) derived in the next section.) 

At the end of this section, we address three technical remarks. First, we discuss the 
validity of the assumption mentioned above (fl8l) . We do not have a mathematical proof 
for the claim that the integration path can be shifted as such described there. However, 
this is plausible because we can prove the result (124p in the following manner: It is 
sufficient to consider the case where A(x) is given by A^{x) = x(x — 1) ■ ■ ■ (x — k + 1) with 
arbitrary positive integer k. The corresponding polynomial A(x) becomes Ak{x) = x k . 
By taking the k-th derivative of the identity 



= e n ° s (27) 

^ n! 

n=0 

with respect to s, and setting s — 1, we obtain 

°° -n 

J2 -±—n{n - 1) ■ ■ • (n - k + 1) = n k , (28) 
' n! 

n=0 

which means 

(A k (n(0))) = A k (n ) . (29) 
Next, following the Ito formula, we find from the Langevin equation fl2U|) for <fi 

= a (kA k + k(k - 1) A-i)^ - ^ (kA k+1 + k{k - 1) A k )^ . (30) 

On the other hand, by using the master equation (CQ), we can confirm that the differential 
equation for d (A k (n(t))) /dt has the structure identical with that of the right-hand side 
of fl30|) . This observation combined with (1291) leads to (A k (n(t))) = (A k (4>(t))) . for any 
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k at any t. Since arbitrary polynomials can be expressed as a linear combination of 
flUI has been proved. 
Second, when n does not satisfy the condition ( {TBI , the integration formula ( |T9l) 
is not available. Therefore, the result (124"|) with the Langevin equation (1231 is valid only 
for the case ({TBI . 

Finally, let us consider two-time quantities such as (y4(n(£i))A(n(£ 2 )))- In order 
to calculate such a quantity, we need a path-integral expression of the conditional 
probability P(n, r|no,0). If we replace the Poisson initial condition with 5(n — no), 
the path-integral expression in this section becomes complicated. 

4. Density fluctuations 

In order to have a direct correspondence with density fluctuations in the coherent-state 
path-integral expression, a non-linear transformation in the action ( FlTl) 

= z/e _x , i(p = e x (31) 

has been employed [8,12-14], which is called the Cole-Hopf transformation. (In the 
standard formulation with (z,z*), (<p,i(f) in (1BT1) is simply replaced with (z,z*).) Quite 
formally, substituting (IBTj) into fflTl) . neglecting terms with third and higher powers of 
X, expecting cancellation of several terms appearing in contributions at t = r, and 
avoiding considerations on the integration path, one can obtain the dynamical action of 
a Langevin equation for the variable v/Vl. This Langevin equation coincides with that 
obtained by system-size expansion methods. 

Although the final result is plausible, it seems difficult to judge the validity of the 
procedures. In particular, in the standard formulation with (z,z*), there is no complex 
number x that would satisfy z = ue~ x and z* = e x , as easily checked for an example 
z — 1 + i. In the formulation with (<j),i(p), the integration path of x is described by 
X = log \ip\ + i7rsgn(y?)/2 with — oo < if < and < <p < oo. The calculation after that 
seems complicated. Toward a justification of ( 1311) . recently, an operator version of the 
Cole-Hopf transformation has been presented [15]; but its mathematical foundation is 
not obvious. 

Based on these understandings, we propose a framework in which the calculation 
procedures mentioned above are formulated without any mathematical difficulties. Our 
basic idea is to introduce a new decomposition of unity associated with the Cole-Hopf 
transformation. As a preparation, motivated by (l3Tj) . we define 

z = z/\z\ (32) 

for any non-zero complex number z. In the argument below, the complex variable z is 
always connected to two real variables u and 9 as 

2 = /ie" i8 , (33) 

where u > and — 7r < 9 < n. Then, z = e~ l9 . By using coherent states parametrized 
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i r°° 

y- 



.in 

oo 



d/i 



e-y £ e- i{£ - m)e \e)(m\ 



2tt 



z) (z\e 



(34) 



We shall employ this decomposition of unity for constructing a path-integral expression. 

We consider the transition probability P{nf, r|ni, 0), which is the probability of 
finding rif particles at time t = r provided that there are n- x particles at t — 0. This is 
expressed as 

1 



P(n f ,r|ni,0) = — {n { 
rif. 



-Ht 



(35) 



Its path-integral expression on the basis of (1341) is written as 

( T POO PIT \ 

n/ d ^/, <" f|2 - > 



t=At 



z t \e 



-HAt 



zt-At)e 



zolm) 



(36) 



(See Appendix A for the corresponding path-integral expression in terms of particle 
numbers instead of coherent states.) As did in the previous section, we have 



\ z t\ e ' 



-HAt\ 



Zt-At) = (zt\zt-At) exp(-F(z t *, z t -At)Ai), where H is expressed as 



IH-M e 



2 e -2i8 t _ At 



(37) 



We then notice 



t=At 



-Ht+Vt-At+m-&t(8t-8 t -At)+0(At 2 ) 



t=At 

T 

IT 

t=At 



t=At 



-ie t (iJ.t-fj.t-At)+o(At 2 ) 



(38) 



where we have used the estimation t — #t_A — 0(At). Furthermore, we have 
(rii\z T ) = z™ f = ^^fe~ mf9T and {z Q \n\) = e ld ° ni . Substitution of these results into 
yields 



( T POO PIT 

n / d ^ / 



2tt 



- i6» (MO - n i ) + ie T (Mr -«f ) 




P(n { ,T\n b O) = limj || / dp t / ^ ) i-n™ f e- nf e - 5({i ^ } ' { ^ }) 
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x -l^e"^ exp[-<S({i0 t }, {//,})] , (39) 

with 

T 

<S({i0,}, {^}) = \ H &> z t-At)At + i0 t (^ - /i t _ At )] . (40) 

t=At 

Here, we perform the integrations with respect to 6 and 6 T . Noting that the limit 
At — > is taken in the final expression, we obtain 

r°° r ■ ,\n.\ i , 

• (41) 

po=n i ,[i T =n i 

This expression is exact and if one simply replaces i9 t and p t with \t and z/ t respectively, 
the resulting action S({xt}, {^t}) is equal to the one obtained by a formal procedure with 
the Cole-Hopf transformation (1311) . Therefore, we claim that our argument provides a 
mathematical foundation for the Cole-Hopf transformation. 

Now, recalling that Q is a dimensionless volume of the system, we focus on the 
regime rif ^> 1 and n- x 3> 1 under the assumption Q 3> 1. More explicitly, by setting 
Pi = n-JVL and pf = n { /Q, we assume a large deviation property 

P(pA r\ Pi n, 0) ~ e -^ (ft ' T|ft ' 0) , (42) 

where p x and pf are the particle densities at initial and final times and are finite in general. 
Note that the relation A(Cl) ~ B(tt) in this paper means (log A(Cl) — log B(Cl))/Cl — > in 
the limit — > oo. The quantity J-(pi, r\pi, 0), which is called a large deviation function, 
characterizes fluctuation properties of a density in a macroscopic system. The problem 
we consider here is to derive a simpler stochastic system that reproduces !F(pf , r\ P i, 0) 
defined in (ji2|). 

Since /io and /i r are fixed to p;f2 and pff2 respectively, the dominant contribution 
of the p t integration in fT4Tl) comes from a region p t ~ O(fi). On the other hand, since 
there is the term exp(i9 t (pt — Pt-At)) in the integrand, the dominant contribution of the 
8t integration comes from a region 9 t ~ 0(1/Q). From these observations, we make the 
transformation of variables as p t = Pt/Q and 7r t = Q6 t in the path- integral expression 
( 14T]) so that the integration over a region where pt — 0(1) and TT t ~ 0(1) provides the 
dominant contribution. We also note the relation ~ n™ f e _nf valid for ^> 1. We 
then have 

-A* ^oo foe Air \ 

(43) 




~ ii m TT / dpt ^ e -s(W,{M) 



t=At 

where S ({iir t } , {p t }) is expanded in terms of 1/fl and At as 

1 



P0=Pl,PT=Pf 



S({m}, { Pt }) = (At) J2 

t=At 

with 



m t ^ 1 (p t , Pt-At) + ^(vr 4 ) 2 ^ 2 (pt-Ai) + R(m, Pt-At) 



*i(pt, Pt-At) = Pt A< - {apt- At - Pp 2 t -At) , (45) 



(44) 



At 



Two Langevin equations in the Doi-Peliti formalism 



11 



*2(Pt-At) = Cipt-M + Ppt-At , (46) 

R{in,Pt-M)=o(^)+0(At). (47) 

The n t integration in f|4*3|) with R(vfft, Pt-At) m <~> being ignored yields the formula 

HPU r\pu 0) = lim min (At) £ ^i^f , (48) 

where min Pi ^ Pf represents the minimization with respect to {/^} with the boundary 
conditions p = p\ and p T = pf fixed. If one includes the term R(m t , Pt-At) in <S and 
treats it as a perturbation, one can show that it does not contribute to the large deviation 
function T . Therefore, as far as we are concerned with the large deviation property, we 
may neglect R(i7T t , Pt-At) hi S. 

Then, by applying the same procedures as (jT9l) and (12TT) in the previous section to 
the expression ( 1431) and (14^1) . we arrive at the final result 

p(p i n,T\p i n,o)^(5(p(T)- Pl )) p , (49) 

where ( ) represents the average over trajectories given by the Langevin equation 



Jt P = *p-Pp +]J—^-V (50) 

under the initial condition p(0) = p\. This result coincides with the result obtained by 
any system-size expansion methods for the master equation. Since (}4"9~|) holds for any r, 
the variable p defined by p/Q corresponds to the density in this description. 

5. Effects of diffusion 

In the main part of this paper, we have focused on the case where diffusion processes are 
sufficiently fast. This formulation can be extended to cases with spatially heterogeneous 
fluctuations. The simplest example in such systems is given by a one- dimensional lattice 
model, which we shall explain. 

Let A = — 0, 1, • ■ • , L} be a one-dimensional lattice. The reaction X ^ X + X 
occurs on each site in the lattice, and a particle moves to a nearest neighbor site at 
rate d. The state of the system is specified by a set of particle numbers on each site, 
n = (n\, • • • , til)- Boundary conditions are assumed appropriately, depending on the 
situation we consider. We then write the master equation for P(n,t) and transform it 
to the equation for a single vector. The time-evolution operator H in the equation for 
this vector is given by 

H = d (dj — a])(dj — dj) + [— a(aj — l)a|dj — (3(1 — a]) a] a? , (51) 

{ij)GB ieA 

where B represents a set of all the nearest neighbor pairs, and dj and dj are introduced 
on each site in a manner similar to d^ and d in section [2j 

Now, as a straightforward extension of the formulation discussed in section HI we 
define a real variable Pi(t) which can be identified with the particle number on a site i 
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at time t. Our central question for this model is to derive an effective stochastic system 
which is defined as the simplest description reproducing large-distance and long-time 
behavior. More explicitly, the effective stochastic system might be defined by a Langevin 
equation which provides a large deviation functional for the transition probability of a 
density field. Since the large deviation functional is related to an effective action in 
the field-theoretical language [19,20], the question is equivalent to an identification of 
the fixed point in renormalization group flow. However, as far as we know, a concrete 
calculation based on such a formulation has not been reported yet. 

Putting aside developing such a theoretical framework, we here present a 
conventional derivation by considering a special situation. We first assume that one- 
site corresponds to a coarse-grained cell in which diffusion processes are sufficiently 
fast so that chemical components can be thought to be uniform in the cell. Let Ax 
be a dimensionless physical size of the cell by setting a reaction length as unity (it 
thus corresponds to Q in the previous sections). As discussed in the first paragraph 
of section [2} we then replace (3 in (1511) with j3/Ax under the assumption Ax ^> 1. 
We also assume that there exists a density field p(x,t), < x < LAx, that satisfies 
p{iAx,t) = p,i(t)/(Ax) and sup^ \d*p(x, t)\(Ax) k = 0(e k ) with a small positive constant 
e. The latter condition means that the density field is slowly varying in space. Since the 
longest wavelength of fluctuations is O(LAx), e corresponds to 1/L. We thus consider 
the case L ^> 1. By repeating the argument in section H] with the replacement of Q with 
Ax, one may derive 



0)^/n*(p(<Ax,r)-^) 



P(n fin , T|n""\ 0) ~ ( | | S [ p{iAx, r) - ^ ) ) , (52) 

where ~ represents the equality valid asymptotically in the limit Ax ^> 1, and ( ) 
represents the average over trajectories of a space-discretized form of the spatially 
extended Langevin equation 

dp 



0f ocp - Pp 2 + dd 2 x p + ^ap + Pp 2 ■£ 1 + d x (y 2dp ■ 6 J , (53) 

under the initial condition p(iAx, 0) = n- mt /Ax. Here, we defined d = d(Ax) 2 and £i 
and £ 2 as noises satisfying 

t)^(x', 0) = 5 lJ 5(x - x')5(t - t') . (54) 

We note that the discretized form of (1541) is 



(£i(£Ax, t)^(£'Ax, t')) = 8 id 8 eA ,-±-8(t - 1') . (55) 

Although these arguments are standard, the validity of the derivation is restricted. 
Formally speaking, we focus on the limiting case L — > oo, d ^ 0, j3 —>■ 0, d//3 2 ~ 
a ~ 0(1) for the lattice model (jSTj) . It is stimulating to analyze (j5ip together with a 
renormalization group idea so as to derive a hydrodynamic description for general cases. 
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Appendix A. Path probability of particle number 



The probability of a particle- number trajectory {n(t)} is also expressed within the Doi- 
Peliti formalism. We briefly explain the expression in this Appendix because there are 
no references mentioning the result explicitly. (See [14] for a related discussion.) We 
start with an expression for the probability P n (t): 



/ T — At oo 



'r-At 



lim 

At^O 



II E \\K{n t ^n t+M ) )P ni {t) 



J=At n t =0/ 



t=0 



(a.i; 



no=ni,n T =n,{ 



where K(n — > m) is equivalent to the previously-defined transition probability (1351) 
applied for a small time step At: 

1 



K{n — > m) 



m 



-HAt 



n 



(A.2) 



and we used the most elementary decomposition of unity 1 = Y^=o \n)(n\/n\. Rewriting 
K{n — > m) as 



K[n — > m) = 



(n + a) 



n + a e 



-HAt 



n) 5, 



tt in v 
aa c -i6(m-n) ^ 



2tt ^— ' (n + crl! 

a=— oo v ' 



m,n+a 



e Wa (n + a 



-HAt 



n 



we express the transition probability during a small time interval At as 



Kin — > m) 



w d9 



-i9(m -n)-^2 L ( ie > n > a ) At 



where L(i9,n;o~) was defined as 

i 

L{\6, n; a) 



n + a 



H 



n > e 



\6a 



(A.3) 



(A.4) 



{n + a)\ 

The expression flA.ll) with flA.31) and flA.41) provides the probability of a trajectory 
{n(t)}. On the basis of this expression, one may develop a system-size expansion as we 
did in section HI This formulation might be regarded as a sophisticated version of the 
discussion in Ref. [2]. 
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